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Abstract 



, First examples of quasi-exactly solvable models describing spin-orbital interaction are 

V/"") ' constructed. In contrast with other examples of matrix quasi-exactly solvable models 

discussed in the literature up to now, our models admit infinite (but still incomplete) sets 
of exact (algebraic) solutions. The hamiltonians of these models are hermitian operators 
of the form H = -A 2 + Vi(r) + (s • \)V 2 (r) + (s + 1) • hV 3 (r) where Vi(r), V 2 (r) and V 3 (r) 
OO ' are scalar functions, 1 is a vector of the angular momentum operator, s is a matrix-valued 

ON . vector spin-operator and h is an external (constant) vector magnetic field. 



1 Introduction 

cu : 

Quasi-exactly solvable (QES) problems are distinguished by the fact that only some of their 
energy levels and corresponding wavefunctions admit explicit construction. In the last decade 
a big progress has been acieved in elaborating the concepts of the phenomenon of quasi-exact 
solvability and formulating methods of constructing and solving QES models of a variety of 



types (for more detail see e.g. the reviews Q, |J, |5)] and book [|§]). In this paper 
we undertake a new step in this direction and construct new classes of QES models which 
(up to now) have never been discussed in the QES-literature. These are QES models with a 
spin-orbital interaction. The potentials of such models have the following general form 

V(x,y,z) = Ui(r) + (s • l)U 2 (r) + (j • h)V 3 (r) 



where r = \J x 2 + y 2 + z 2 is a radial coordinate, 1 = (l x ,ly,lz) = (i(yd z — zd y ),i(zd x — 
xd z ),i(xd y — yd x )) is hermitian (vector) operator of angular momentum, s = (l x , l y , l z ) is her- 
mitian (vector) spin-operator whose components realize a certain unitary finite-dimensional 
representation of algebra so(3), j = s + 1 is the operator of total momentum, and h = 
(h x , h y , h z ) is an external (constant) magnetic field interacting with the total momentum. 

One of the most important distinguishing features of QES models which we intend to 
present here is that they have an infinite number of exactly (algebraically) constructable 
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energy levels and corresponding eigenvalues. They are however quasi-exactly solvable because 
the set of their exact solutions is still incomplete and does not fill all the spectrum of a model. 
First examples of such models were presented in our recent work [Q] where they have been 
called "infinite QES models ". 

Another unusual feature of these models is that they (in contrast with models discussed 
in paper [Q]) are matrix models with physically realistic hermitian hamiltonians. Everybody 
who has some experience with quasi-exact solvability in the matrix (multi-channel) case knows 
how difficult is to satisfy the condition of hermiticity when constructing such models. It is 
hardly neccessary to remind the reader that up to now only a couple of hermitian matrix 
QES models have been constructed (see e.g. [§, [§]]). 



2 Starting point 

To demonstrate how does our construction procedure work we start with the simplest one- 
dimensional QES model with hamiltonian 

H = + (c ~ 1/2 ^ C ~ 3/2) + [b 2 - 2a(2m + c + l)]r 2 + 2abr 4 + aV (1) 

acting in Hilbert space of functions defined on the positive half axis r € [0, oo] and vanishing 
sufficiently fast at its ends r = and r = oo. Here a, b, c are real parameters satisfying the 
conditions (a > 0, c > 0) and m is a non-negative integer. As it was demonstrated in [||]], 
for any fixed m the Schroedinger equation 

Hif>(r) = Eip(r) (2) 

for model [l] admits algebraic solutions whose general form is given by the formulas 

tf(r) = r c ^ II (r 2 /2 - &) exp (-^ - \\ (3) 

m 

E = 26(2m + c) +8a^& (4) 
i=l 

The m complex numbers £j in expressions ^ and || satisfy the system of m algebraic equations 

E 7T7 + ^- 6 -H=0, i = l,...,m. (5) 

k=l,k& & Sk 

It turns out that system || has only m + 1 permutationally invariant solutions for any given m 
which are represented by the sets of real points £j. Each solution is completely characterized 
by a (quantum) number k = 0,1, ...,m which indicates the number of positive ^-points. 
According to formula ||, the number of positive ^-points determines the number of (real) 
wavefunction zeros, which, in turn, determines the ordinal number of an excitation (oscillation 
theorem) . This means that model [l] has m + 1 exactly constructable solutions describing the 
ground state and m first excited states. A more detailed exposition of properties of model [l] 
and its algebraic solutions can be found in the book [M]. 
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3 The modified equation 



It is not difficult to see that the transformation 

ip(r) = r<p(r) (6) 

reduces the equation [2] to the form 

+ U m) ) ^ = Ev{r) (7) 

in which we used the notation 

V(r, I, m) = [b 2 - 2a(2m + 5/2 + l)]r 2 + 2abr 4 + a 2 r 6 (8) 

and 

Z = e-3/2 (9) 

Hereafter we shall consider I as a new independent parameter taking (by agreement) only 
non-negative integer values. The form of the first three terms in the equation coincides 
with the form of the radial part of a tree-dimensional Laplace operator. For this reason it 
seems quite natural to interpret I as the 3-dimensional angular momentum and try to relate 
the equation ^ to a certain 3-dimensional quantum problem. In the following three sections 
we show that there are three such possibilities leading to three different kinds of quasi-exactly 
solvable problems in the 3-dimensional space. 



4 The first possibility 

One of the simplest possibilities of interpreting equation [7| is based on the assumption that 
function ^ entering into [7| is I independent: 

V{r, I, m) = V (r, N) = [b 2 - 2a(N + 5/2)]r 2 + 2abr A + a 2 r 6 (10) 

For this the number 

N = l + 2m (11) 

must be fixed. In this case, equation [?] takes the form of a typical radial Schroedinger equation 
for a spherically symmetric 3-dimensional equation 

(-A + V (r, AO) y, z) = EV{x, y, z). (12) 

Since both m and I are assumed to be positive, the condition [ll] leads to a finite number of 
possibilities with m = 0, 1, [N/2], and I = N,N — 2, N - 2[N/2], respectively. For this 
reason, for any given N, the model |l^ is quasi-exactly solvable and has (as usually) only a 
finite ( [N/2]([N/2] + l)/2) number of explicit solutions. The model of such a form and even 
its more complicated spherically non-symmetric versions were considered many years ago in 
papers [©,§]. 
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5 The second possibility 



Another possibility of interpreting equation p]is to consider m as a fixed number not restricting 
the value of I. In this case the function || becomes linearly dependent on I and can be 
represented in the form 

V(r,l,m) = Vi(r,m) - I ■ V 2 (r) = {[b 2 - 2a(2m + 5/2)]r 2 + 2abr A + a 2 r 6 } - / • {2ar 2 } (13) 



It is quite obvious that, in order to associate the equation 13 with a certain 3-dimensional 
Schroedinger equation, we must find a proper 3-dimensional source for the term which is 
linear in the momentum I. The first think which comes in ones head is to look for the 3- 
dimensional scalar operators O which wuld commute with both the Laplace operator and r 
and would have the eigenvalues linear in I. In this case we could consider as a reduction of 
a 3-dimensional problem 

(-A + Fi(r,m) - O ■ V 2 (r)) V(x,y,z) = EV(x,y,z) (14) 

The linearity in I means that the operator must be proportional to the operator of the angular 
momentum 

1 = (l x , l y , l z ) = (i(yd z - zd y ),i(zd x - xd z ),i(xd y - yd x )). (15) 
But this is a 3-dimensional vector while the operator we are looking for must be a scalar. 



The only possibility to construct a scalar from 15 is to take a scalar product of 1 with another 



vector operator. It is quite obvious that there is no such operator if we restrict ourselves to the 
single-channel problems. However, if we admit the consideration of multi-channel problems, 
then a good candidate for the second operator can immediately be found. This is obviously 
the spin operator s! Rerstricting ourselves (for the sake of simplicity) to the 1/2-spin case (2 
by 2 matrices), we can easily check that the spectrum of the operator 

<3 = 2-sl (16) 

(which, obviously commutes with both A and r) is linear in I. Indeed, representing operator 
|j~6| in the form 

0=j 2 -l 2 -s 2 (17) 

(where j = 1 + s is a total momentum) and taking for concreteness a particular case with 
j = I + 1/2 we easily find the corresponding branch of the spectrum 

o = j(j + 1) - 1(1 + 1) - s(s + 1) = (I + 1/2)(Z + 3/2) - 1(1 + 1) - 3/4 = I. (18) 

This finally leads us to a 3-dimensional matrix QES models 

(-A + {[b 2 - 2a(2m + 5/2)]r 2 + 2abr 4 + a 2 r 6 } - 2(s • 1) • {2ar 2 }) f(x, y, z) = E^(x, y, z) 

(19) 

describing spin-orbital interaction. 

It is a time to ask ourselves of what kind of models did we obtain? First of all, one should 
stress again that these models are actually quasi-exactly solvable. This follows from the fact 
that for any given m and I they have an infinite number of normalizable solutions, but only 
m + 1 of them are exactly (algebraically) constructable. Second, and this is may be the most 
important thing, despite the fact that the set of exactly constructable solutions is incomplete, 
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this set is infinitely large. This is so because the number I is not fixed by the 3-dimensional 



model |19| It appears as a solution of the eigenvalue problem for operators O and may take 
arbitrary non-negative integer values. 

In conclusion of this section note that the hamiltonians of models we obtained are hermi- 
tian by construction. 



6 The third possibility 

The last interesting possibility of reducing the equation ^ to a 3-dimensional form appears 
when the function || depends on both parameters I and m. In this case the function |8| becomes 
linearly dependent on both I and m and can be represented in the form 

V(r,l,m) = Vi{r)-l-V 2 (r)-m-V 3 (r) = {(b 2 -5a)r 2 +2abr 4 +a 2 r 6 }-l-{2ar 2 }-m-{4ar 2 } (20) 

By analogy with the previous section we can consider the numbers I as the eigenvalues of the 
operator of spin-orbital interaction, and the only thing which remains to do is to interpret 
m as an independent quantum number appearing in equation [7| as an eigenvalue of a certain 
operator M commuting with the variable r, Laplasian A and the spin-orbital operator s • 1. 
A good candidate for such an operator is the z-projection of the total momentum s + 1. In 
fact, it should not neccessarily be a z-projection. Because of the spherical symmetry, it could 
be equally weel a x- or y- projection, or any other projection. We can therefore represent 
this operator in a covariant form 

M = 2(s + l)-h (21) 

where h is a unit magnetic field. We introduced an additional factor 2 to make the eigenvalues 
of operator M integer rather than half integer. Of course, the negative integers are not 
interesting for us, because, as we remember, only for non-negative integer values of m the 
system admits algebraic solution. Summarizing, we can consider [7| as a reduction of a 3- 
dimensional problem 

(-A + Tq(r) - (s • 1) ■ V 2 (r) - 2(s + 1) • h-V 3 (r)) *(x, y, z) = EV(x, y, z) (22) 

which can be treated as a spectral problem for a matrix quantum model describing spin-orbital 
interaction together with the interaction of a total momentum with an external magnetic field. 
It is remarkable, that the model ^ does not contain any integer parameters anylonger. All 
these parameters appear dynamically as solutions of the eigenvalue problems for additionally 
introduced symmetry operators. At the same time, the model ^ remains quasi-exactly 
solvable, because for any particular values of these eigenvalues the equation has only a 
certain incomplete set of solutions. 



7 Conclusion 

The method of construction infinite (matrix) QES models exposed in this paper is, obviously, 
quite general and can easily be used for building other spin-orbital models with more com- 
plicated potentials and higher matrix dimensions. For this it is sufficient to start with other 
known one-dimensional QES models first rewritting them in the form of a radial Schroedinger 
equation and then interpreting the /-dependent terms appearing in their potential as the 
eigenvalues of a spin-orbital operators. 
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